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We review some motivation behind the introduction of chiral random 
matrix models in QCD, with particular emphasis on the importance of the 
Gell-Mann-Oakes (GOR) relation for these arguments. We show why the 
microscopic limit is universal in power counting, and present arguments 
for why the macroscopic limit is generic for a class of problems that defy 
power counting, examples being the strong CP and U(l) problems. Some 
new results are discussed in light of recent lattice simulations. 



1. Introduction 

It is well known that the spontaneous breaking of chiral symmetry in 
QCD along with its explicit breaking, provide stringent constraints on on- 
shell amplitudes whether at threshold or beyond. These constraints can 
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be exploited in power counting (ChPT) Q, or to any order ||, leading to 
a comprehensive description of a number of reaction processes. The same 
symmetry constrains the off-shell amplitudes such as the vacuum to vacuum 
amplitude in the presence of external sources. 

What is less well-known perhaps, is the fact that in a finite Euclidean 
volume V = L 4 similar constraints are still in action even in the case where 
the box size L is smaller than the pion Compton wavelength (in 
units where the QCD scale A = 1). This is a regime where the spontaneous 
breaking of chiral symmetry is obsolete. This observation for the off-shell 
amplitudes was first emphasized by Gasser and Leutwyler Their ob- 
servation carries to most models with spontaneously broken symmetries in 
d > 2 as was discussed by Hasenfratz and Leutwyler || . 

The regime for which mV > 1 will be referred to as macroscopic. In 
it, chiral symmetry is spontaneously broken with (q'q) = 1 (again in units 
where A = 1). The regime for which mV < 1 will be referred to as micro- 
scopic. In it, chiral symmetry is effectively restored with (q'q) = mV + O. 
The nature of the corrections O will be discussed below. The regime mV ~ 1 
will be referred as transitional. In it, chiral symmetry is about to be re- 
stored or broken depending on how we wish to look at it. Throughout, the 
box size is chosen L > 1/A = 1, so that nonperturbative effects are present. 
At small but finite temperature, a similar classification holds in Euclidean 
space with V = f3L 3 . 





Fig. 1. QCD in a symmetric (left) and asymmetric (right) box. 



The importance of the parameter mV for the constant mode problem 
in QCD was first noted by Jolicoeur and Morel § in the strong coupling 
regime. This remark is important as it implies that the observations to 
follow carry for both strong and weak coupling provided that a continuous 
symmetry is broken spontaneously with the occurrence of Goldstone bosons. 

In the present work, we will go over a number of observations in the 
context of chiral symmetry that paves the road for the onset of chiral random 
matrix models in QCD with their relevance to the three regimes of mV 
detailed above. In section 2, we go over the original observation by Banks 
and Casher fn. In section 3, we detail some of the arguments given by Gasser 
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and Leutwyler Q regarding the finite volume partition function, paying 
particular attention to the Gell-Mann-Oakes-Renner (GOR) relation. In 
section 4, we summarize how these observations were used by Leutwyler and 
Smilga []|] to derive spectral sum rules, and derive new spectral sum rules 
for generalized GOR relation. In section 5, we present simple arguments 
for why these sum rules are minimally reproduced by chiral random matrix 
models. In section 6, we go over the construction by Verbaarschot and 
Zahed || for the microscopic n-point correlation functions as the master 
formulas for all diagonal and off-diagonal sum rules. The n-point correlation 
functions specify the degree of flavor mixing in the n-vacua, and uniquely 
condition the bulk spectral rigidity and level variance. The relevance of these 
observations for lattice simulations in symmetric and asymmetric boxes is 
discussed. In section 7, we discuss the microscopic spectral density in the 
double scaling limit, following the original arguments by Jurkiewicz, Nowak 
and Zahed |1C] and more recently Damgaard and Nishigaki [11| and Wilke, 
Guhr and Wettig 111. New sum rules with decoupled sea and valence quark 
effects are discussed. We argue that the associated spectral density in the 
double scaling limit is universal in QCD, within a range of quark masses 
that we discuss. In section 8, we put forward arguments for why chiral 
random matrix models are of interest in the macroscopic regime, especially 
when chiral power counting breaks down. We illustrate our points for the 
CP and U(l) problems. In section 9 we summarize our conclusions. 



2. Quark Spectrum 

In the early 80's Banks and Casher made a remarkable observation re- 
garding the relation of the quark condensate to the QCD quark spectrum. In 
an Euclidean volume V, they noticed that the amount of quark pair conden- 
sation (q'q) relates in a simple way to the distribution of quark eigenvalues 
z^(A). Indeed, 

J A + im 



w)) 



-in v\m) 



where the limit is understood for mV > 1. The quark eigenvalues de- 
pends implicitly on the given gauge configuration, through the eigenvalue 
equation iD q^ = XkQk- They are paired, since [y^,iD] + = 0. This pairing 
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Fig. 2. Quark spectrum in the vacuum. 



causes z^(A) to be an even function on its support. The averaging in ([l]) is 
over the gauge configurations. 

In (|l|) the zero virtuality point A = [] is the analogue of a Fermi surface. 
This observation suggests in fact that the regime mV > 1, in which chiral 
symmetry is spontaneously broken, is actually diffusive [14|, following the 
derealization of the quark eigenmodes. 

More remarkably, however, is the fact that for the pair condensate to 
be non-zero, the density of quark eigenmodes near zero virtuality has to 
grow with the volume V. In other words, the level spacing of the quark 
eigenmodes in a finite volume V and near zero has to scale as 1/V in sharp 
contrast to the level spacing of a free particle in the same box which has to 
scale as 1/L. 

The spontaneous breaking of chiral symmetry causes a large accumula- 
tion of quark eigenmodes near zero virtuality as pictured in Fig. for a fixed 
gauge-configuration. The occurrence of exact zero modes of the nonzero n 
(winding number) configuration is geometrical and follows from the topo- 
logical character of the background gauge configuration. Conversely, any 
model that breaks spontaneously the underlying chiral symmetry displays a 
spectrum of the type shown in Fig. g. This means QCD, instanton models, 
NJL models, etc. 



3. Finite Volume Partition Function 



3.1. With GOR 



Few years later, Gasser and Leutwyler Q noticed that in a finite volume 
V for which the Compton wavelength of the pion is large, that is mV < 1 
but L > 1, while chiral symmetry is effectively restored with (q'q) = mV + 
O, the way it is restored is still uniquely conditioned by the way chiral 



1 The quark virtuality A is analogous to a complex mass in 4d Euclidean space, or an 
energy in l+4d Minkowski space. 
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symmetry is about to be spontaneously broken. For QCD this observation 
is natural if we were to approach the microscopic regime mV < 1 from the 
macroscopic regime mV > 1 by fine tuning the bare quark mass m. 

mV > 1 region: In a finite Euclidean box V the finite volume partition 
function Z is the sum of zero-point motions of singlet excitations: pions, 
kaons, nucleons, etc. (Casimir effect). In the regime mV > 1 the spon- 
taneous breaking of chiral symmetry together with the Gell-Mann-Oakes- 
Renner (GOR) relation [15] allows for a book-keeping of InZ by using chiral 
power counting with a typical momentum scale p, where U 1/L ~ p 

and m ~ m 2 ~ p 2 and m?V ~ 1 fixed. In this book-keeping, the contribu- 
tion to order p 2 is given by the following Lagrangian [|| 



F 2 I 
C 2 = — Tr (d^d^U) - -S Tr(mU + £/W) 



4 V p p ) 2 

for Nf flavors with U = e llT l F (tt = vr°A Q ) an SU(Nf) valued field subject 
to the boundary condition U(x + L) = U{x). The normalization is chosen 
TrA^A^ = 25 Qj/ g. Throughout we will use m for either the non-degenerate 
case m = (mi, m^f ) or the degenerate case m = mljv r We hope that 
the reference to each case will be clear from the notation used. 
In a periodic box, the free pion propagator is 

g ih n x ~ g ih n x 

G -( X ) =yY, k 2 +m 2= y^2 +yY, k 2 + ^2 ( 3 ) 



where the primed sum is over non-zero four momenta with k n = (m, n2, n.3, 11.4) 
in units of 2tt/L. To order p 2 the pion mass is m 2 = mT,/F 2 ~ m and the 
condensate is —(qq) = S. To order p 4 , the one- loop contribution from (J2j) 
has to be supplemented with the most general terms consistent with broken 
chiral symmetry and general principles Q (Lorentz invariance is upset in a 
box to order p 4 ). 



\ Si v n 

\ It t \ 
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/ \ » 1 I 



Fig. 3. Meson contribution to InZ for mV > 1. 



mV^ < 1 region: In the microscopic regime 1 > mV, the zero mode 
contribution 1/mV in (||) dwarfs the non-zero mode contributions (primed 
sum) causing the chiral power counting to break down. In this regime Gasser 
and Leutwyler suggested a reorganization of the chiral power counting fixing 
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mV instead of m 2 V with expansion in e = 1/L, thereby allowing for a sys- 
tematic resummation of the pionic zero modes in the finite volume partition 
function. 

To order 0(e 2 ) the dependence of Z[S,P,V] on the external sources 
S = mS /2 + S' and P is simply 

Z[S, P, V] = N[0] [ dU e ^Tr((5+ 4 P)t/+h.c.) _ (4) 
JSU(Nf) 

N[0] receives contribution from the non-zero mode terms, essentially the 
Casimir energy of free pions in a box V as indicated by the first diagram 
in Fig. 3. This contribution is S, P-independent to the order quoted. For 
x = rnYiV and Nf = 1, Z[m] ~ e x , while for Nf = 2 with equal masses 
it is Z[m] ~ I\{2x)/2x 0]. For one massive flavor and Nf — 1 massless 
flavors it is Z[m] ~ lN f -i(x)/x N f^ 1 . In each case, Z[m] = Z[m,0,V] and 
the proportionality constant is m independent. 

From (fD the chiral condensate to order 0(e 2 ) is 

itfq) = [ dC/Tr(C/ + C/t )e ^VT r([ / +f 7t) (5) 

Z[m\ 2 Jsu(N f ) 

and vanishes as x = mV — * 0, to the exception of Nf = 1 where i(q*q) = 
— S. For Nf = 2 and two equal masses, i(qiq) = -S/ 2 (2x)//i(2x), while 
for one massive and Nf — 1 massless quarks i(q^q) = —al^ f [x)/lj<i f -x{x). 
In this limit all points on SU(Nf) are equally weighted. Although chiral 
symmetry is restored in the regime mV < 1, the way by which it is restored 
is conditioned by : 

i) the nature of the coset; 

ii) the explicit (Nf,Nf) breaking; 

iii) the underlying GOR assumption. 



+ X # X + XXXX + 



Fig. 4. Zero momentum contribution to hiZ for mV < 1. 



3.2. Without GOR 
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mV > 1 region: The GOR assumption is important in establishing 
@) as the starting point for power counting in ChPT in the macroscopic 
regime. This assumption may be relaxed in the form of generalized ChPT. 
An example to order p 2 is 

F 2 1 
£ 2 = — Tr [dprfdpU) - -VE Tr(mU + C/W) 

+ATi{m)Um ] U + h.c.) + B ( r Er(m t i7 + h.c.)) 2 

+C (Trirn^U - h.c.)) 2 + DTx{m)m) (6) 

as suggested by Stern, Sazdjian and Fuchs and others fl6|| . Here E and 
m are counted of order p, while A,B,C,D are counted of order p°. The 
constants E, A, C are low-energy parameters, while D is not (see below). 
To order p 2 the pion mass squared and the chiral condensate are 

m 2 w = ™ (E - 8m^ - 8mN f B) 
- (qq) = E - m(2A + L> + AN f B) . (7) 

The constant D is subtraction dependent. It relates to the m-dependent 
quadratic divergence in the condensate. It is not amenable to low-energy 
constraints. 

mV < 1 region: Since the contribution of the pionic zero modes in (|3|) 
are of the form 1/Vm^ 2 , the finite volume partition function (Q) for fixed 
Vm 2 would involve A,B,C,D from @. The microscopic limit is changed 
compared to the GOR assumption if we were to take Vm 2 fixed but small 
(see section IVB). Such a limit could be reached by varying m for fixed 
m 2 V < 1. The power counting in e = 1/L would be: E and m of order 
e 2 , and A,B,C,D of order e . In this regime, the chiral condensate again 
vanishes in the chiral limit as all points on the SU(Nf) manifold are equally 
weighted. 

Clearly mV fixed selects the GOR form in the microscopic regime, and 
could be used for a precise determination of E by comparison to ChPT. This 
method may be more accurate than a simple extrapolation of the condensate 
to zero m, owing to finite size effects. This is more than an academic 
exercise since both (||) and (JgJ) lead to the same low energy constraints, 
and experiments are underway at DAPHNE to test the validity of the GOR 
hypothesis by firming up the accuracy of the S-wave tttt scattering length 
near threshold | |17| |. 



4. Microscopic Sum Rules 
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4.1. With GOR 

In the early 90's Leutwyler and Smilga ||] remarked that the finite vol- 
ume partition function could be used to derive sum rules for the quark 
spectrum (qualitatively shown in Fig. ^) in the microscopic regime for gauge 
configurations of fixed winding number n. Assuming the concept of n-vacua 
for the QCD state, the finite volume partition function for fixed n relates 
to the finite volume partition function through 

Zn [ m ] = / _ e~ ine Z[me ie / N f] 
Jo 2ir 

= N(0) [ dU (detU) n e !VETt(mtwW) (g) 

Ju(N f ) 

in the microscopic regime mV < 1. For QCD n = ^ f E a ■ B a and is valued 
in Z. The left-hand side in (|8|) involves the quark spectrum averaged over 
gauge configurations with fixed winding number n, while the right-hand 
side involves an information from the hadrons. Z n [m] is local in the quark 
variables, but non-local in the hadronic variables. Scalar n-point functions 
following from (||) do not obey the cluster decomposition. 

Z n [m] does not discriminate between Nf = 1 and Nf > 1 in QCD. Also, 
Z n [m] is insensitive to the nature of the sources S' and P in (|j) , and only 
depends on the quark mass matrix chosen diagonal, and its overall phase 
for n / 0, since || 

Z n [VmW] = (det(VW))"Z n [m] . (9) 

This is not the case for (0). 
For one flavor, we have 

Z n [m] =m |n| ( J] (Xl+m 2 )) n , A =n\l n (x). (10) 

A fc >0 

with x = mVY, and I n a Bessel function (the primed product does not 
contain the exact zero modes and the eigenvalues with negative real part). 
The averaging is over the gauge configuration A with fixed winding number 
n. Using the infinite representation for the Bessel function and rearranging 
(|10D it follows that 

«?K)».-9( i+ £) 
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where £ nj fc are the zeros of the Bessel function J n (£, n ,k) = with k = 1, 2, .... 
The double averaging on the left-hand side of ( |IT| ) now includes m' n l FJl -^1 
as part of the measure, with a normalization to 1 at n = and m = 0. By 
identifying the same powers of m on both sides of ( [HI) sum rules for the 
quark eigenvalues at fixed n follow, in the microscopic regime. An example 
is | 

where we have again reinstated £ and Nj for convenience. The contribution 
of the ultraviolet modes to the primed sum in (12) vanishes like V/V 2 in 
the large volume limit. 

It can be checked that the sum rule (^) is just a suitable arrangement of 
moments directly amenable to the finite volume partition function (Q) . They 
reflect therefore on the sum rules satisfied by the invariant QCD correlators 
at zero momentum in the microscopic limit. For example 

1 ' 1 Y 2 r 27T riff / ie ~ iB \ 2 

^«E4»» = ^i £**{(W* + Ue"f)))„ (13) 

where the averaging on the rhs is done with the Haar measure normalized 
to 1. Notice that while the lhs is local in the quark variables, it is not on 
the rhs, owing to the 9 integration. The rhs is the zero momentum and 9- 
averaged contribution from the quark susceptibility (9\(q'q) 2 \9), evaluated 
in the microscopic limit mV — * 0, with (q'q) = 0. The susceptibility reduces 
to the variance of the scalar operator on the invariant measure. 

The remarkable property of Qi~TD is that it implies that the quark spec- 
trum near zero virtuality is indeed 1/V spaced, with a spacing that follows 
as if from a 'master' gauge field in the form £ n k/V. In the microscopic 
regime the continuum limit V — > oo is taken by keeping the level spacing 
finite. The master formula for all diagonal microscopic sum rules is just 
given by 

-. OO £ oo 

v n {s = \V) = -J2S(X-^) = J2 6(8 - Cn, fc ) • (14) 

k=l k=l 

The limit V — > oo with s fixed was originally discussed by Shuryak and 
Verbaarschot j^], and the microscopic sum rules were numerically checked 
to hold for the instanton liquid model. This is expected, since condition 
i)-iii) at the end of section IIIA hold verbatim for the instanton liquid 



model as shown by Alkofer, Nowak, Verbaarschot and Zahed [19j. 
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4.2. Without GOR 

We note that the present results were derived starting from the GOR 
assumption with the e expansion carried in the regime mV fixed but small. 
If the GOR relation is given up as discussed in section III, then in the regime 
m^V ~ m 2 V fixed the microscopic sum rules are changed. In particular ( |l2|) 
becomes for n = 



l»fln_ fV 2 1 (\TrU\\ D 

The averaging on the rhs is done over U(Nf) with the Haar measure. The 
1/V term is of order e 4 , but so is X 2 . We note, that contrary to m which 
is a parameter, £ is fixed by the QCD dynamics. In small volumes V with 
m^V held fixed but small, the power counting is valid, and the second term 
may not be negligible. The D dependence in (|i~5| ) is expected from the 
subtraction dependence to this order. 



5. Chiral random matrix models 



With the renewed interest in the instanton liquid model along the lines 
discussed by Diakonov and Petrov 



and also Shuryak | 21[| , it became 
clear that the bulk issues elegantly revealed by the variational and statis- 
tical approach to the instanton physics are generic and likely amenable to 
a description in terms of random matrix theory [22]. This observation is 
particularly relevant for the quark eigenvalues and their correlations as orig- 
inally suggested by Nowak, Verbaarschot and Zahed p3| ]. In this context, 
the relevance of the three Wigner ensembles p5] was noted earlier by Si- 



monov [24], and the relation to the microscopic limit made later by Shuryak 
and Verbaarschot 18 ] . In this limit the physics is that of a rotor plus a po- 
tential for the uncolored Goldstone constant modes Q , and that of random 
matrix theory for the colored constant modes. 



5.1. With GOR 



The universality of the finite volume partition function in the micro- 
scopic regime implies that any model that satisfies the conditions i) — iii) 
summarized at the end of section IIIA, will do as far as the microscopic 
sum rules are concerned, thanks to the observation made by Gasser and 
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Leutwyler [|j . The minimal model that satisfy these requirements is a four- 
fermi model in O-dimension (matrix model). Specifically, 

C = 9 t imeWNf q + J^q R q L q R q L + ... (16) 

with d\m.qi R = Nj x til,r are Grassmann numbers, 75 = (l R , 1l), and 
n L + n R = V . The ... terms in (|l6|) could be either singlets (1,1), or non- 
singlets. The latters should vanish in the chiral limit as m a with a > 1. 
The finite volume partition function associated to ( |T6| ) may be either 

cast in the Grassmann variables, or by bosonization through A ab = qL a Qjt 
which is NfriL x NfUR, turned into an integral over bosonic matrices A in 
the form 

Z[6,m}=[ dAe- y ^ A ' A ^ + - 

ime i0/N f A \ , , 

The determinant is V x V for Nf = 1. From our preceding arguments, it 
does not matter what ... is in the microscopic limit mV < 1. The finite 
volume partition function (|17|) is that of a chiral random matrix model. In 
QCD the quarks are in the complex representation, so A is a complex valued 
matrix. The chiral random matrix ensemble generated by the description 
(0) is ChGUE. 




5.2. Without GOR 



The minimal four-fermi model in O-dimension that is commensurate with 
©is 



,t 



vme 



iy 5 9/N f 



q R qLq R qL 



+ 



(B + C) 



V 



\A[q\rn)q L f + (q R mq R )' 



+2 



( \q R qR)(qL mm ^qL) 



v 

+DTv{rm)m) + ... 



(18) 



Again, ... are either singlets or non-singlets but then of higher order in m. 
The finite volume partition function can be bosonized using: A = qiq R , 
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B R = QrQr an d B L = QlQli the result is a multi-random matrix model 
(N f = 1) 



ime ie + me i9 B R A 

ime~ id + me~ id B L 



with width of the distribution given by 



where 



(19) 



B L ab B L cd ) = (BX) = J= • ^^ls ad 6 bc (20) 



V 

A abA\ d ) = ^ • —^—5 ad 5 bc (21) 

1 V JEFF 



^ = ±J 1 + ( B-cf-£) (22) 

JEFF ^ V ^ / 



The inverse bosonization of (18) using q\nqR,L as bosonic variables can be 
shown to yield (0) as expected. In this case, the different scalings in V in 



the width (21) are important. The terms of order (TrC/) 2 follows only after 
an exact gaussian integration over A and B has been performed. We note 
that in the present formulation, (|l8|) yields A in terms of B, C and X. 



6. Microscopic Spectral Distribution 

Since the microscopic sum rules follow from chiral random matrix models 
in the microscopic regime mV < 1, we can then use the powerful machinery 
of random matrix theory to make precise statements on the eigenvalue dis- 
tribution and its correlation. Without loss of generality, we can just assume 
a gaussian weight in (fl7|) and proceed. Throughout, we will discuss the 
GOR case. The alternative case is straightforward to implement. 



6.1. Orthogonal Polynomial Method 

Following Verbaarschot and Zahed Q, we may specialize to the symmet- 
ric case ul = ur = V/2 = N (zero winding sector), and use the polar de- 
composition of the unitary matrix A, A = UAV*, where A = diag(Ai, Ay) 
in R+, and U in U(V) and V in U(V)/U(l) v . The division by U(1) N pre- 
serves the number of degrees of freedom. As a result, the angular and radial 
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integration separates, the former inducing only an overall m-independent 
normalization. It will be ignored. With this in mind, we have 

POD JV J 

Z[N, m}= J] d\ie V ^ Y[(X 2 + ™/) A W (23) 
Jo i=i / 

with the Vandermonde determinant 

A[A] = niA,|(A?-A J 2 ) 2 (24) 

i<j 

following from the Haar measure on SU(V) in the polar decomposition 
[p5| . The first determinant in (|2^) stems from the integration over the 
Grassmannians. 



In random matrix theory [22|, the partition function ( J23| ) has a simple 
form in terms of the joint eigenvalue distribution 

i/(Ai,...,A v ) = |(Pi...Pv|Ai...Ay)| 2 (25) 

where the P^s form a set of ortho- normalized polynomials of degree % — 1 
with 1 < i < V (with the term of degree i normalized to 1), with the 
measure 

/ dX 2 e~ vx J] (A 2 + mJ)Pi(A)Pj(A) = l id . (26) 
Jo /=1 

The relation ( |25| ) follows simply from the properties of determinants (in- 
variance under arbitrary linear combinations of either columns or rows). 
It is a Slater determinant with Pi,...Py as the 'bra's and Aj, ...Ay as the 
coordinate 'ket's. In particular, 

r v 
u(X) = / dX 2 ...dX N v(X, A 2 , ...Ay) = IP(A)| 2 (27) 

i=l 

for the eigenvalue distribution at finite V and iVy\ Eq. (25) summarizes all 
eigenvalue correlations. 

In the microscopic regime m < 1/V. Consider the simplest case where 
m = 0, for which (p6| ) is just the condition for the generalized Laguerre- 
Selin polynomials |J7 In the microscopic regime V — ► oo with VA fixed, 
this distribution can be evaluated exactly in the form 

u„(x) = {T?x/2) (j 2 Nf {Y,x) - J A r /+ i(Ex)J Ar/ _i(Ex)) (28) 
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in the sector with zero winding number. The dependence of v s (x) is shown 
in Fig. ||, for Nf = 0, 1, 2. The larger Nf the wider the distribution at zero 
virtuality, following larger repulsion from the fermion determinant. In our 
case £ = 1, but in general it is the scale associated to (q'q) in the chiral 
limit as is apparent from (||). 

In a finite box, (q'q) gets renormalized by non-zero modes, and thereby 
develops a dependence on 1/L along with m (and of course S) as originally 
discussed by Hansen and Leutwyler ^6|. We will refer to it as E(m,L). 
Its dependence on the low energy parameters at the chiral point such as 
the pion decay constant, scattering length etc. along with 1/L and can be 
organized using a double book-keeping in e and p. 

1.6 i ■ ■ ■ 1 



> 



1.2 



0.8 



0.4 



I . . . 1 

5 10 

x 

Fig. 5. Microscopic spectral distribution for ChGUE with Nf = (solid), Nf = 1 
(dashed), Nf —2 (dotted). 

The correlations between eigenvalues in the microscopic limit can be 
constructed using similar arguments. In particular, the connected correlator 
between two eigenvalues x and x' in the microscopic regime for n = is |§ 

v a {x,x') = Y?xx' (29) 

(xJN f (T,x)JN f -\(Y l x l ) — x' J]\f f (T,x') Jiv / _i(Ex) 
x 2 — x' 2 

The connected two-point correlator (|2^ ) obeys the consistency condition 

J dx v s (x,x') = (30) 

which reflects on the conservation of the total number of eigenvalues. We 
will see below that the class of all n-point unconnected correlators may be 
resummed to give the microscopic eigenvalue distributions with m / in 
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the double scaling limit. Eq. ( p9| ) can be used to check the off-diagonal sum 
rules discussed by Leutwyler and Smilga ||. The microscopic two-point 
correlator (|29|) may be used to address some issues related to the flavor 
admixture and spectral rigidity in the microscopic limit, as we now discuss. 



1 









10 



Fig. 6. Microscopic connected correlator for ChGUE Nf = (solid) and Nf 
(dashed line) with an arbitrary normalization. 



As we have stressed above, the universality of fl28|) irrespective of the 
gaussian measure adopted is guaranteed by i) — iii) as discussed in section III, 
thanks to the modified power counting. These observations are of course in 
agreement with those made by Nishigaki and Damgaard |27| in the context 
of polynomial weights, though more general, as they apply to a field theory. 
Although the present arguments apply to QCD in even-dimensions they can 
be easily modified to accommodate QCD in odd dimensions as discussed by 
Verbaarschot and Zahed [28|. This point is noteworthy as it implies that 
the issue of the instantons and a vacuum angle is not crucial for the present 
observations. 



6.2. Flavor Mixing 



It was noted early on by Nowak, Verbaarschot and Zahed [23|, that the 
eigenvalue correlator in the macroscopic limit reflects on important aspects 
of the QCD vacuum as a disordered medium. In particular, it was suggested 
using the instanton liquid model that the admixture of flavor (strangeness) 
is suppressed in the ground state (average of all n-states) if the instanton 
density is closer to the metallic regime. A measure of this suppression is 
given by the two-point correlator in the eigenvalue distribution, and could 
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be directly related to the level variance in the macroscopic limit. Zweig's 
rule is supported by Wigner-Dyson statistics in the macroscopic limit. 

In light of our present discussion, similar questions may be asked in the 
microscopic limit, for which it is easy to show that the Zweig's violating 
correlator investigated in |23|| , that is (s^su^u) can be written in terms 
of (|2^ ) in the n = state with yi = rriiV as, 

i 2 (shv)u) (y u ,y s ) = dxdx' Vu Vs v s (x,x'). (31) 
J x*+y u z x u +y s A 

In particular, 

i 2 (s^ sv) u) n X 4 , . 

lim — s ifl = 32 

y-0 y u y s lW 2 {N f + l) V > 

which is seen to break the cluster decomposition. The admixture of flavor 
in the microscopic limit in an n-state is universally given by the multipoint 
correlators in the microscopic limit. 



6.3. Spectral Rigidity 

The microscopic spectral density u s (x,x') characterizes the level spacing 
distribution a quantum leap apart, not only at zero virtuality but also in 
the bulk of the spectrum. Indeed, it carries important information on the 
spectral statistics in the microscopic regime, such as the spacing distribution 
or the level variance. It also plays a major role at the interface of the theory 
of disordered systems and the semiclassical theory of quantum chaos, as we 
will try to elaborate further on in the end. 

To exemplify some of these points, let us consider in more detail the 
genealogy of S2, the variance in the number of single particle levels within 
an energy interval A A centered around Ao, in the microscopic regime Ax = 
VA(X — Ao) 3> 1 but fixed as V — > 00. This regime is characterized by 
constant mean density, v s {x) = ^/V,^] and in the decomposition of the full 
two point correlator, u s (xi,X2) = v s {x{)v s [x2)[8{xi — X2) — R(xi,X2)], the 
cluster function R(x\,X2) is given as 



R(s,A) 



AsinSs + (-iy v /scos 2SA 

AS7T 



(33) 



2 Actually this assumption is not really required, as unfolding procedures of spectra 
may be used |p2|. 
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with A = {x\ + X2) /2 and s = x\ — x 2 - In the bulk (Ao 7^ 0) the second term 
may be neglected and the cluster function depends only on s |29|]. The edge 
(Ao = 0) is special due to the vanishing of the density of states as required by 
chiral symmetry. The level number variance, T,2(N) is expressed through 
the cluster function R as 

[■N 

E 2 (N)= dxxdx 2 R(s,A). (34) 
Jo 

Specifically 

Z 2 (N) = -^logN (35) 

with (3 = 2 in the bulk and twiced at the edge, Ao = 0. For the orthogonal 
and symplectic case (3 = 1 and 4 respectively. Both for the quenched and 
unquenched cases, the spectral rigidity is 21nA r /(/?7r 2 ) as opposed to N for 
Poisson statistics. The decrease in the variance is an indication of a more 



rigid spectrum due to level repulsion 3C]. 



6.4- Spacing Distribution 



The probability P(s) to find two energy levels in the quark spectrum a 
distance s apart in the microscopic limit, that is s = (Aj — Xj)V/A fixed 
as V — > 00, is also intimately related to the properties of the two-point 
correlator (^). For s <C 1, that is as Aj — * \j, we immediately see from the 
Vandermonde determinant (|24|) that P(s) ~ s@ with (3 = 2. This is the usual 
level repulsion between two neighboring levels as predicted by Wigner [|30| . 
It follows from the random lore for 2x2 matrices, since in this limit the 
rest of the eigenvalue spectrum decouples. Clearly it is independent of the 
number of flavors, except when one of the level is at zero virtuality. In this 
case, further repulsion is introduced by the fermion determinant and the 
winding number n in the massless case as is evident from Fig. H 

For s S> 1 (still in the microscopic limit), the spacing distribution is 
solely governed by the asymptotics of the two-point correlator (p9|). For 
R(s) ~ l/(3s 2 , the level distribution asymptotes — lnP(s) ~ 0s as was 
shown by Dyson [29] using the the Coulomb gas treatment for (Q) in the 
quenched case. In the unquenched case, the asymptotics of the two-point 
correlator is unaffected by the fermion determinant, resulting into the same 
level distribution. Indeed, for any Nf, it is straightforward to see that 
the fermion determinant contributes to the potential term (one-body) not 
the interaction term (two-body). Since it is the latter that conditions the 
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level spacing distribution, the self-quenching of this result is immediate. 
Hence, the whole level spacing distribution for (§) could be described by 



the standard Wigner surmise [3C] 



P(s) = As p e~ B ^ (36) 

where (3 = 2 and A, B are pure numbers. Which is drastically different from 
the Poisson distribution P{s) = e~ s , for uncorrelated states. 



6.5. Relevance to Lattice 



By now the present observations have been generalized to the case of 
finite winding number, as well as the two other Wigner ensembles: ChGOE 
and ChGSpE by Verbaarschot 31 ] , and checked against lattice Dirac spectra 
by Schafer, Weidenmiiller and Wettig |52], and others |33|]. In regard to the 
lattice measurements, we would like to point at few issues in regard to the 
strong coupling aspects of the simulation as well as the approach to the weak 
coupling regime. Some of our remarks will also relate to certain aspects of 
the quenched approximation. 

In strong coupling, the Kogut-Susskind action is characterized by an 
exact U{1) x U(l) symmetry. The analogue chiral random matrix model is 
ChGSpe [Ml. So the strong coupling regime is characterized by universal 
microscopic spectral oscillations, even though Lorentz and full chiral sym- 
metry are still absent. We also expect quantitative changes in the sum rules 
for quenched and unquenched simulations since the character of the coset 
and hence the number of Goldstone mode changes in this case ||. 

In weak coupling, the Kogut-Susskind action is expected to undergo a 
transition to a phase with continuum Lorentz and full symmetry. A struc- 
tural change from ChGSpe to ChGUE is therefore expected, for otherwise 
the true continuum would not have been reached. This point can be checked 
at the level of the microscopic sum rules, with due attention to scale trans- 
lations, e.g. V,m (see below). 

In symmetric boxes an accurate measurement of the finite volume chiral 
condensate in weak coupling can be reached by checking the microscopic 
sum rules, or simply monitoring the Bessel oscillations. This means an 
accurate assessment of the chiral condensate at order p 2 , as well as the pion 
decay constant to the same order through possible finite volume effects. 

In asymmetric boxes, the chiral condensate gets modified by tempera- 
ture effects which can be assessed using a double expansion in e and p (see 
below) at low temperature (small asymmetry). At high temperature, the 
chiral condensates disappears as the level spacing near zero becomes larger 
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than 1/V. For temperatures T ~ T c , universality arguments may be used 
to argue for critical sum rules [Q]. For T > 3T C , the virtual spectrum is 
gapped by the lowest Matsubara mode A ~ irT, except for thermodynami- 
cally irrelevant zero modes. 

Finally, the relevance of the present concepts to lattice QCD in the 
strong coupling indicates that they are also applicable to lattice electronic 
systems with ferromagnetic or antiferromagnetic ground-state order. An 
example is the undopped Lanthanum-Copper-Oxide material for high tem- 



perature superconductors [35]. 



7. Double Scaling Regime 

7.1. The Formulae 

In the transition regime mV ~ 1 the current quark masses affect the 
distribution of quark eigenvalues quantitatively for Nf ^ 0, through the 
occurrence of the fermion determinant. Surprisingly enough, the non-local 
character of the latter supersedes the effects of local weights to upset the 
conventional universality arguments based on polynomial weights in the 
microscopic limit. This point was originally shown by Jurkiewicz, Nowak 



and Zahed [10] using the supersymmetric method for Nf = 1 and zero 



winding number, and more recently by Damgaard and Nishigaki [11] and 
Wilke, Guhr and Wettig [jl^] using the orthogonal polynomial method with 
arbitrary Nf and winding number n 3 . 

The massive spectral density in the double scaling limit yields by defini- 
tion the quark condensate in a fixed winding number configuration n, but 
with different m s sea and m v valence quark masses in general. Specifically, 

i(q*q)n(£,y)= dx u sn (x,y) = y) (37) 

Jo £ + 4 

with the rescaled valence quark mass £ = Em„7 and a finite sea quark 
mass y = T,m s V. We are using different notations for £, y to highlight their 
different origins. For £ = y the result is just the massive quark condensate 
in a winding number configuration n. The physical quark condensate with 
^ = y follows by suitably averaging over all windings n. For Nf = 1 and 
n = the massive spectral density is (x = T,VX) 

l^u sfi (x,y)= f[.#(*) + J?(*)] (38) 



3 Incidentally, this construction applies to any additive part one adds to the fermion 
determinant say T provided that it is suitably rescaled T —> VT fixed, in the micro- 
scopic regime. 
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J (x) 



[yJo(x)h{y) + x,Ji(x)I (y)] 



I (y) x 2 + y 2 

for x > 0. For y = @ reduces to (||) with 2V/ = 1. 




2468 02468 10 

x x 
Fig. 7. Massive microscopical spectral density from an ensemble of 50000 N=50 
random matrices (dotted line) for m s = 1 (left) and m s = 0.1 (right), Nf — 1, and 
the analytical formula ( |38"1) (solid line). 



7.2. Sum Rules 



The microscopic spectral density (^) yields new microscopic sum rules 
with sea quark effects included. These sum rules are not a priori amenable 
to the finite volume partition function discussed above simply because ( |38l) 
selects solely the sea quark mass effect, which is not a physical concept. 
Finite volume partition functions for QCD with different sea and valence 
quark masses can of course be simply constructed by using quark ghost 
fields ||. Such discriminations are easily achievable on the lattice, and 
hence are of theoretical relevance. An example being the detailed valence 
quark study undertaken recently by Chandrasekharan and Christ [|F|, for 
a fixed sea quark mass contribution. 

The derivative of (|37| ) with respect to the valence quark mass generates 
a string of sum rules that are now sensitive to the effects of the sea quark 
mass y. Specifically, for n = 0, Nf = 1, and fixed £,y 

^«£^»o = (39) 

where the double averaging on the left-hand side involves m' n ' + 
as part of the measure, in contrast to the one used in (|l2| ) with m s = 0. 



z98ii printed on February 1, 2008 



21 



£o(£, V) depends on the rescaled valence and sea quark masses indepen- 
dently. Its analytical form is 



£ 2 



+ 



2^o(OL w yh(y) 



(40) 



y 2 -e 



MO 



For y = it is in agreement with the result discussed by Verbaarschot 
In the limit of equal sea and valence quark masses the expression simplifies 
to 



E 2 



(41) 



In Fig. |a we show the behavior of So(^,y)/S for y = (plus signs), 
y = (, (crosses), y = 5£ (triangles) and for fixed m s = 1, that is y = 1/N 
(boxes) from numerical simulation of 5000 N = 50 random matrices and the 
corresponding analytical form ( |40| ) (lines). For fixed and finite sea quark 
mass Eq. ((3^) diverges for small £ as — log £. In Fig. |8]b we show the behavior 
of the n = condensate from (pg|) for the same set of masses, namely y = 
(solid line), y = £, (dashed line), y = 5£ (dotted line) and fixed m s = 1 
(dashed-dotted line) . 

1 



A 
cr 
cr 

V 




0.5 




10" 



10' 







10 z 




Fig. 8. a Sum rule (39) obtained numerically from an ensemble of 5000 N = 50 
random matrices and the analytical result, b The quark condensate for various 
masses. See text. 



For two massive sea quarks the result for the microscopic spectral density 



was worked out by Damgaard and Nishigaki [11], and Wilke, Guhr and 
Wettig for arbitrary n. In this case, the analogue of (pQ) is 



£ u (g,2/i,2/2) 
£ 2 



I K + I X K X 



(42) 
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-vl 



{ 



yiifi b -y2is4yt(yl-e)+yl(e-y 2 i)H 4 (y!-y 2 2) 

y 1 y 2 I?I b 1 [(yl-y%)I K +(yl-e)lSKg+(e-y 2 i)I b K b 

M. yi qi b + y2 i b i$) [ay!-y 2 2)ioK 1+yi (y 2 2 -eKK? 

e 2 ( y i 2 - y2 2 )/i^+ y i 2 (2/ 2 2 -c 2 )/i a ^+ y2 2 (e 2 -^)/W]} 



with Ii = If = I{(yi) and l\ = h{y2)i respectively and similarly for 
K^s. For one massive quark and Nj — 1 massless quarks, the sum rule is 
given by 



E 2 



(43) 



y 



+ 



y 2 I Nj+ i(y) 



I Nf -iK Nf+ i - I Nf -i(y)K Nf+ i(y) 



N f I N x {y){y*-?) 



iNf-lKjNf-V 



-I Nr -i(y)K Nf -i(y) 



For one massive quark and Nf — 1 massless quarks, the sum rule for the 
massive quark is S 2 (y,y) = ^ 2 In s /2yJjV/— l for n = as also noted by 
Damgaard [O]. We note that in this case Z = Zq, which is the partition 
function in the n = sector. In particular, i{q i q) = i(q^q)o- Additional 
sum rules involving diagonal and off-diagonal eigenvalue correlations are of 
course possible. The off-diagonal ones will involve the n-point correlation 
functions of eigenvalues with massive quarks. 
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Fig. 9. Normalized condensate. See text. 

Recently Chandrasekharan and Christ 12] have analyzed in details the 
behavior of the valence quark condensate for fixed sea quark mass m s a = 
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0.01, for staggered Nf = 2 QCD, over several decades of the valence quark 
mass in an asymmetric lattice V = 4 x 16 3 a 4 . Their analysis was carried 
for various lattice coupling (3 = 6/g 2 (varying lattice spacing a), around 
the chiral transition point (3 C = 5.275. Their underlying quark spectrum is 
clearly sensitive to temperature, and an analysis of this point using a macro- 



scopic model distribution was discussed by Nowak, Papp and Zahed [p7 |. 
For small asymmetries, (3 < f3 c , we may still assume that the conditions 
i)-iii) summarized at the end of section IIIA are still valid, in which case 
i(q^q)o{Ci y) may be amenable to a lattice comparison. For that, we need 
to identify E^a 3 for different j3 and the proper dimensionless combination 

= toZj^V on the lattice. For the former, we have: S^a 3 = 0.2217 
(J3 = 5.245), E L a 3 = 0.1357 (J3 = 5.265) and E L a 3 = 0.0543 (J3 = 5.270). 

In the left part of Fig. |9| we show the lattice results for the normalized 
condensate versus the rescaled valence quark mass, The line is the zero 



flavor result [36] in the n = topological sector seemingly in good agreement 
with the two flavor lattice data (symbols for different /3's). In the right we 
show the two flavor result ( [42"D for equal valence masses, y = m„EK. The 
solid lines represent Nf = 0, 1, 2 massless flavors |3(|, respectively, the plus 
signs are for the lattice data with (3 = 5.245 while the crosses for (3 = 5.270. 
We conclude that the used sea mass on the lattice is still too large making 
the simulation analogous to a quenched one (Nf = 0). Decreasing the sea 
mass by a factor of ten we start to see the deviation from the quenched 
result (boxes), and lowering the mass by a factor of hundred we recover the 
Nf = 2 zero mass limit (triangles). The same effect can be achieved dialing 
(3 closer to the chiral transition point, where the condensate, E, disappears. 
This tendency is seen in the right part Fig. ^ for [3 = 5.270. 

For small valence quark mass £ < 0.1 limit of Eq. (|40|), we find 

In the large sea mass limit we recover the quenched limit , — log £/2 — C + l/2 
with a logarithmic divergence. 

7.3. Universality 

The result (^) is universal for the QCD spectrum as we now show. 
Following Smilga and Stern [38], we can rewrite the massive fermion deter- 
minant for a fixed background as (n = 0) 



2 

m 



A m = II ( A I + m ") Nf = A oexp ( N f £ ln(l + ^-) ] . (45) 



A fe >0 \ A 



k 
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Hence, the density of eigenvalues for finite m reads 

1 / roo m 1 \ 

p(X,m) = -(explNfJ^dX' u(\',A) ln(l + — )J u{X,A)) A . 

where the averaging is over A including Ao, the fermion determinant with 
zero mass quarks, and u(X, A) is the unaveraged spectral operator 

u(X,A) = J26(X-X k [A]). (46) 

k 



In terms of the rescaled variables x = VX and y = Vm, the expression ( |4q 
reads 



u s (x,y) = (exp ^NfJ dx' u(x',A) ln(l + 



which involves all the moments of the microscopic spectral distribution. 
Specifically, 

u s (x,y) = v a (x) + N fJ Q dx' In ^1 + -^^ (lo(x',A) oj(x,A))a + ••• 

which is the weighted density-density correlator (second moment) in the 
microscopic limit. Each of the expectation value is carried in the n=0 state 
with massless quarks. All the moments are universal and given by random 
matrix theory. For Nf = 1 the result for fli^ ) is (|38|). 



7-4- Range of validity 

Since mV ~ 1 is the transition region in which the quark condensate 
(||) is no longer averaging to zero, the pionic zero modes in (|3|) are no 
longer dominant. What is the range of validity of ( |38| ) in QCD? The answer 
follows by noticing that (^) can be extended to one extra dimension, that is 
a Lagrangian in 1+4 dimensions. The support on the fifth direction is [0, /3] 
with periodic boundary condition for U (xs + (3,x) = U (x$, x). Here plays 
the role of a 'temperature' and — kiZ/fl for large f3 is just the ground state 
energy of the singlet Hamiltonian described by 

f*2 TnV 

Hi+4 = W ~ —Tr(U + C/t) (47) 
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where L 2 is the Laplace-Beltrami operator on the SU(Nf) manifold. Note 
that the present construction is analogous to the one discussed by Leutwyler |3_ 
in 1 + 3 dimensions. 

For mV = the spectrum is that of a spherical top with the irreducible 
representations of SU(Nf) as eigenf unctions. The spectrum is 1/V spaced 
with a mass gap. The ground state wavefunction is a constant on SU(Nf). 
The first excited state is Nj degenerate. For mV < 1 the degeneracy is 
lifted. For Nf = 2, the degeneracy is 4 with a triplet (pions) and a singlet 
(sigma) state. With increasing mV, the triplet states are pushed down and 
the singlet state up. The former merge into the physical pion mass as shown 
schematically in Fig. |l0[ We note that m < 1/V puts the mass range into 



the ergodic regime, and 1/V < m < 1/yV to the diffusive regime [14|. 
The transition to the macroscopic regime with a mass gap given by ^fm 
(remember that m\ = 2m with our conventions) sets in for mV ~ 10. 



E„-E n 



1/V 



1/V 



1/1/ 



Fig. 10. Mass gap versus m from (pTT) 



7.5. New Results 



Recently Damgaard [11] has suggested an interesting relation for u s (x, y) 
in terms of the finite volume partition function involving 1 + 2 flavors with 
respective masses (y,ix,ix), using a property of the chiral random matrix 
measure. Generically, the microscopic spectral density reads [11] 



p(X) = C ■ |A|n(A 2 + mi) 

i 

Schematically it means that 



2 Z(mi, 



Z(m 1 ,...,m Nf ) 



(48) 



- A)) = (n(A. - ^) 2 ) • |A| n(A 2 + 



(49) 



i.e. we may trade in a Dirac delta function for two additional flavors with 
imaginary masses. The origin of this formula in the random matrix model 
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is very simple. One has just to insert the Dirac delta function into the ex- 
pression of the RMM partition function written in terms of an integral over 
eigenvalues. Then the two flavors and the factors outside the expectation 
value just follow from contracting with the Vandermonde determinant. The 
term with the probability distribution e ~ NV ( M ) does not contribute in the 
microscopic limit. 

One would like to understand this formula from the point of view of 
QCD. The main difficulty is that in this regime of QCD, we are dealing only 
with the Goldstone boson degrees of freedom and all the explicit dependence 
on the spectrum of the Dirac operator seems to have been lost. In particular 
the low lying states are composite states in terms of the quark fields. 

The main unsolved problem is to find a nonlocal operator in the mesonic 
picture whose expectation value would correspond to the microscopic spec- 
tral density [40]. 



8. Macroscopic Regime 

In the macroscopic regime mV > 1, the constant modes in (p) cease to 
have a preferential role, and we are back to full QCD with (0), (|6) (or even 
other unexplored power countings) as starting points in power counting. 

However, there are a number of problems in QCD where power counting 
(in the sense exposed) breaks down. Examples are the strong CP problem, 
the U(l) problem, and phase transitions to cite afew. Since these problems 
involve in an intricate way the issues of large volumes, small quark masses, 
zero modes etc. it is useful to address in a framework where these quantities 
are simply separated in a way that allows for an analytical treatment. More 
importantly, the framework should be able to include some generic dynam- 
ical aspects of QCD in the form of few vacuum moments and symmetry. 

The chiral random matrix models introduced above for studying the 
microscopic character of Dirac spectra, offer such an example. In many 
ways they can be regarded as a schematic description of the chiral physics 
at work in a cooled lattice QCD configuration within a finite Euclidean 
box V. Since these models allow for closed form results in terms of m, 
V, Nf, etc., they are very useful for addressing some subtle aspects of the 
thermodynamical limit in the presence of zero or near-zero modes whether 
in vacuum or matter. 



8.1. Wide Correlators 



Chiral correlations in QCD involve usually eigenvalues which are sepa- 
rated by a macroscopic distance in the Dirac spectrum. They are defined 
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for V — > oo and fixed Aj — Xj. Correlations over these distances will be 
discussed in this part of the lecture. For bounded spectra, Ambj0rn, Ju- 
rkiewicz and Makeenko [41] have shown that the smoothened n-point corre- 
lation functions could be classified by the support of the spectral densities, 
independently of the specifics of the random ensemble and genera of the 
topological expansion. In QCD the Dirac spectrum is not bounded due to a 
strong ultraviolet tail p(A) ~ |A| 3 . This tail, however, is not very important 
for infrared physics, except for multiplicative renormalization factors and 
anomalies. In a lattice formulation, this tail can be subtracted through a 
cooling procedure. The resulting spectrum is likely bounded, (except for 
possible tails which may be triggered by a partial loss of confinement), in 
which case the results discussed by Ambj0rn, Jurkiewicz and Makeenko EJ 



for the hermitian case, and Janik, Nowak, Papp and Zahed [42] for the non- 
hermitian case may apply. An early account may be found in the analysis 
by Nowak, Verbaarschot and Zahed of the instanton liquid model [23]. 



8.2. U(l) problem 

An important problem in QCD relates to the fact that the rj in na- 
ture is much more massive than the 7r, K, rj system. It is believed that the 
discrepancy in mass is related to the fact that the U(l) current in QCD 
is anomalous. Indeed, using Ward identities it follows that (Minkowski 
space) 



%o P = - jy3<??} + J d x ( T *Vo(x)vo{0)) (50) 
where the topological susceptibility reads 

Xtop = / d 4 x(T*Z(x)Z(0)) (51) 



and the topological density is 8irZ(x) = a s E a ■ B a (x). The singlet current is 
rjo(x) = qi^x,q{x). For small m a gap in the singlet correlator requires that 
Xtop 7^ 0. This is the presently accepted view for the resolution of the U(l) 
problem, although some difficulties may be noted |f44f. 

The chiral random matrix models discussed in (|17|) offers a simple way 
to model some of the aspects of the present problem. Indeed, the matrix 



ime A 
ime~ 



(52) 
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with entries ni and ur can be viewed as a schematic description of the 
topological zero modes in a finite volume V, with J 3 = — ur. Since 75 
is just 



75 



1 
-1 



(53) 



it follows that the form of the anomaly in the present model is just Tr75 = 
nL — ur. The non-vanishing of ( |5l| ) in QCD, will appear naturally in this 
model by summing over matrices of the type ( |52[ ) as in ( |i~7| ) with varying 
sizes. The sizes will be gaussian distributed according to 



e 2v x * 



(54) 



Clearly {{til — nn) 2 ) = V\* with the determinant set to 1 in ( |T7| ) (quenched 
case). 

This model has been discussed recently by Janik, Nowak, Papp and Za- 
hed p5[] . Their results are shown in Fig. 11, where the topological suscep- 
tibility ( |Tlja), the pseudoscalar susceptibility (lib) and the quark conden- 
sate ( |Tl|c) are studied versus m for different values of N = til + tir = V . 
Clearly, the extrapolation to small values of m warrant larger and larger 
sizes ./V (volume V) for the result of Fig. |i~l| b to be meaningful. The solid 
curves are analytical results in the thermodynamical limit. The numerical 
results were obtained by sampling over matrices of different sizes using a 
Gaussian distribution for the matrix elements as in ( |i~7|) and the distribution 
(|54D for the size variations. 



8.3. CP problem 

The non-vanishing of the topological susceptibility means that the vac- 
uum partition function of QCD depends on the value of the vacuum angle 
8. This point is also clear from our earlier arguments as the free energy 
F(0,m,V) = — laZ[0, m, V]/V was 9 dependent by construction. QCD 
with at small 9 implies a free energy shift F(9) — F(0) ~ Xtop# 2 /2, which vi- 
olates T and P. Since strong interactions are known empirically to preserve 
T and P in the vacuum, this causes the strong CP problem. 

A variety of scenarios have been put forward to resolve it ranging from 
axions to confinement [46]. This problem could have been in principle settled 
on the lattice if it were not for the breakdown of conventional Monte-Carlo 



algorithms at finite 9. Recently, Schierholz [47] has addressed this issue in 



the context of the CP n model in two dimensions. His conclusions that the 
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Fig. 11. Normalized topological (upper left), pseudoscalar (upper right), scalar 
(lower left) susceptibilities, and Ward identity (lower right) for Nf = 1. The 
numerical simulations where carried for fixed 27V = n + + n_ = 60, 120, 200 on 
an ensemble of 10 4 , 5000 and 5000 matrices, respectively, and < (n + — ra_) 2 >= 
Nx* = N. The solid line is the analytical result p5[ . 



model exhibits a first order phase transition at finite were challenged by 
Plefka and Samuel [48|. 

Our present framework offers a testing ground for the present problem 
since both the issues of quenched, unquenched, finite size and current quark 
effects can be dealt with explicitly. As was shown by Janik, Nowak, Papp 
and Zahed the model in many ways resemble the effective models dis- 
cussed using effective Lagrangians in the saddle point approximation. The 
advantage, however, is that the present model allows for numerical sim- 
ulations that test for the validity of such an assumption, the subtlety of 
the chiral and thermodynamical limit and the importance of the numerical 
accuracy. 



9. Conclusions 

We have reviewed some of the motivations in the introduction of chi- 
ral random matrix models to QCD problems, emphasizing some universal 
aspects in the microscopic limit as well as some generic aspects in the macro- 
scopic limit. In particular, we have shown that by relaxing the GOR relation 
new microscopic sum rules and random matrix models may be set up in the 
context of QCD. 

We have discussed the double scaling regime and we have presented 
relations for the chiral condensate in a fixed winding number sector, that 
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Fig. 12. ChRMM (quenched) free energy for N=250 (left). The normal precision 
calculation is shown by circles, the finite volume high precision by dashed line and 
the large volume high precision by solid line. Similar results are obtained for the 
free energy in the CP™ model (right) ]47| . 

are sensitive to the sea and valence quark masses independently. 

In the macroscopic regime, we have suggested that ChRMM although 
not universal, exhibit generic aspects of the finite volume QCD problem 
that are useful for addressing currently open problems in QCD, with em- 
barrassing similarities with bulk lattice simulations. ChRMM offer useful 
framework for understanding the interplay between the thermodynamical 
limit, the chiral limit and the precision of numerical algorithms. 

Finally, we would like to conclude by pointing out that the chaotic as- 
pect of the quark eigenvalues near zero virtuality as first revealed by random 
matrix theory from the finite volume partition function, may be established 
directly from QCD and general principles without recourse to power count- 
ing These arguments carry to lattice formulations as well as supersym- 
metric QCD. 
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